arXiv: 1501.0477lvl [cond-mat.mes-hall] 20 Jan 2015 


Drag of ballistic electrons by an ion beam 


V. L. Gurevich and M. I. Muradov 

Ioffe Institute, Russian Academy of Sciences, 194021 Saint Petersburg, Russia 

Drag of electrons of ID ballistic nanowire by a nearby ID beam of ions is consid¬ 
ered. We assume that the ion beam is represented by an ensemble of heavy ions of 
the same velocity V. The ratio of the drag current to primary current carried by 
the ion beam is calculated. 

The drag current appears to be a nonmonotonic function of velocity V, it has 
maxima for V near v n p/2 where n is the number of electron miniband (channel) and 
v n F is the corresponding Fermi velocity. This means that the ion beam drag can be 
applied for ballistic nanostructure spectroscopy. 


PACS numbers: 68.47.Fg,73.20.-r,73.23.-b 


I. FORMULATION OF THE PROBLEM 


Drag as a physical phenomenon in solids can be described as follows. Consider a solid 
with two types of quasiparticles (type 1 and type 2). One creates a flux of the quasiparticles 
of type 2, the so-called driving current. As a result of the interaction between particles a 
current of quasiparticles of type 1, the so-called drag current is excited. An example of such 
phenomenon is the Coulomb drag where a current in conductor creates a current in adjacent 
conductor — see papers by Pogrebinski 1 ' and Priced. 

The purpose of the present paper is to consider a somewhat different situation where the 
driving current is created by real particles outside the conductor (rather than by quasipar¬ 
ticles within it). This would provide a contactless method to generate a drag current (or 
voltage) in a nanostructure. 

Two formulations of the problem are feasible. 

1. The dragging flux consists of heavy ions of almost the same velocity V. 

2. One considers a flux of weakly ionized gas that is in thermal equilibrium having some 
temperature T and hydrodynamical velocity V. 


In the present paper we treat the first possibility, 
beam, i. e. a flux of ions having the same velocity V. 
of velocity V is determined by the accelerating voltage 


MV 2 


= e/TJ. 


In other words, we consider an ion 
For the simplest situation the value 
QJ and the ion mass M as 


( 1 ) 


Here ej is the charge of an ion. 

For the drag system we will treat a ballistic (collisionless) electron transport in a quantum 
wire. Such nanoscale systems may have rather low electron densities that can be varied by 



2 


means of the gate voltage. The collisionless quantum wires act as waveguides for the electron 
de Broglie waves. For a strong Fermi degeneracy 


T < p 


( 2 ) 


where T is the temperature (we will use the energy units for it) while /i is the Fermi energy, 
each miniband of transverse quantisation (channel) gives the following contribution to the 
conductance 


Gn = 


7T H' 


( 3 ) 


(e being the electron charge) so that the total conductance is 

G = AfG 0 . 


Here A f is the number of such active channels, i. e. the minibands with bottoms e n (0) below 
the Fermi level /j. 

Our purpose is to investigate the main features of this drag phenomenon. We assume 
that the distance d between the ion beam and the wire is much larger than the width of 
the wire, so that on the scale of this width the Coulomb interaction of ions and electrons is 
a smooth function. Then the selection rules for the corresponding matrix elements require 
that electrons involved in the transitions change their quasimomenta but remain within the 
initial transverse quantized channel n. One can vary the velocity V of the ions with the 
accelerating voltage 03 and measure the resulting variation of the drag current (or drag 
voltage). We will denote by V r the volume occupied by the nanowire while Vr will be the 
volume where the flux of ions propagates and interacts with the electrons of the nanowire. 
We assume both V r and Vr to have a ID shape of length L parallel to z-axis. 

One can give the following qualitative considerations concerning the drag by an ion beam. 
Due to the conservation of such quantities as the energy, the transverse quantized channel 
number n and the (quasi)momentum in the electron-ion collisions one has to consider in 
the Born approximation the transition of electron from | n,p) to \n,p + q z ) state (where p is 
the ^-component of electron quasimomentum) and that of the ion from |P) to |P — q) state 
according to relation 

p!_ = .(p+&) 2 . ( p -q) 2 

2m 


2 M 2 m 2 M 

The ^-function describing the energy conservation can be therefore written as 


( 4 ) 


<5 


r 2 2 

-^-(1 + m/M) + — (p — mV) + 

2m v m v 2 M 


2 m 


I Qz 


5 [q~ - 2 (mV - p )], 


( 5 ) 


where P z = P = MV. Further on we will take into account that m/M <C 1 and ne¬ 
glect m/Ad as compared to 1 and (■ m/M)q\ as compared to q Therefore the transferred 
(quasi)momentum is q z = 2/mV — p) and the probability of such a transition includes the 
factor 


fnp{ 1 fn,p+q z ) fn,p+q z \ 1 fnp) fnp fn,2mV—p 


( 6 ) 
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as well as the electron-ion Coulomb interaction matrix element squared. For the ID situation 
under consideration it has a factor proportional to 


KH\q z \d/h) 


( 7 ) 


q z =2(mV—p) 


where d is the distance between the ion beam and the wire and Kq is the McDonald function 
[see below Eq. (JTFJ)]. One can use for it the following approximate equations 


2 


Kq(s) ps In 


( 8 ) 


for s«l, 



( 9 ) 


for s> 1 


where In 7 = 0.5772. The drag current is proportional to the sum over electron momenta p 
of the products in Eqs. ([5]) — (J7J) . We consider p < 0 and require the state p to be occupied, 
this condition leads to — pp < p < 0. The requirement that the hnal state with momentum 
2 mV — p should be empty gives 2 mV — p>pp provided V < vp /2 (see Fig. [I]) . If V > vp/2 
there is no additional restriction except — pp < p < 0 , i.e. all occupied states are involved 
in transitions. Therefore, if V < Vp/2 we get for the drag current 



2mV-p 


-P p 2mV-P F 


FIG. 1: Momenta from — pp to 2 mV — pf are involved in transitions. If V > vf/2 all negative 
momenta p contribute to the drag current. 



( 10 ) 


and we see that increase in V decreases the minimal transferred momentum and increases 
the effective Coulomb interaction K 0 (2pd/h). If V > vp/2 we have 



( 11 ) 


and increase of V results in a decrease of the drag current. These equations provide an 
adequate description of the drag current dependence on the ion beam velocity as can be 
readily seen in our quantitative approach below. 
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II. INTERACTION OF ION BEAM WITH ELECTRONS OF 

NANOSTRUCTURE 


For simplicity, we assume the width of the beam to be constant (actually it may slightly 
vary in the course of beam propagation). Then one can write the distribution of the ions 
within the beam as 


F P = N(2icH) 3 S(P x )S(P v )S(P, - P), 


( 12 ) 


N being the ion concentration. 

The collision term of the Boltzmann equation for ID electrons and 3D ions is in the Born 
approximation 


df 


pn 


dt 


coll 




V K d 3 P f V R d 3 g 2vr 
( 27 t h) 3 h 


\(p,n,P\U\p + q z ,n,P -q>|' 


^ d(y€ n p + Ep £n,p+q z h/p_q) [/ n p(l fn,p+q z )Fp fnp+q z {1 fnp)Fp— q] (13) 


where 

e n (p) = e n (0) + p 2 /2m 


(14) 


Here n is the number of the channel, i.e. of the mini-band of ID transverse quantization 
(according to the assumption made above this number does not change in the course of 
electron transitions), q is the transferred (quasi)momentum, M is the mass of an ion, m is 
the effective mass of conduction electron and 


U = 


eej 


(15) 


1 + k |R — r| 

describes the Coulomb interaction of the ion with charge e/ and electron in the wire, hi being 
the dielectric susceptibility of the wire. For the matrix element in Eq. (fT3l) we have 

/J '' 2 eej 


(p, n, P\U\p + q z ,n, P - q) = / d A r / d A Ri>* n (r ± )T 

Jv r Jv R 


P L(1 + /t)|r — R| 


Since 


where 


d.Zdz 


e^’-z)/* = 2K 0 (\q g \\Ar ± \/h) 


L\r — R| 

| Arj_| = [x - X) 2 + {y- Y) 2 


i(ri)'f P - q e , ^ /fi 

(16) 

(17) 


we can write 


(p,n,P|U|p + ^,?r,P-q) = 

(i t^Vr / dR± j dr ^( r ^\ 2 e ~ l ^ n±/hK o(M Ar ±\/ h )- 

The Boltzmann equation for electrons is 

Sfnp _ 
v dz 



(18) 

( 19 ) 
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where 


de 


v = 


pn _ P 

m 


dp 


is the electron velocity. 

To calculate the current in the wire, we iterate the Boltzmann equation for the electrons 
of the wire in the term describing collisions between electrons of the wire and ions. In the 
zeroth approximation one can choose the electron distribution function in the collision term 
to be equilibrium one. In what follows f np = fp(e np ~ p) will be implied, where fp is the 
Fermi distribution function while // is the Fermi level. The first iteration of Eq. (1191) gives 
for the nonequilibrium part of the distribution function 


A fnp — ( 2 2 1 v I{f > F} 


( 20 ) 


J np 


for p > 0 and p < 0, respectively. Here I{f,F} is a shorthand notation for the collision 
term. Using the particle conservation property of the scattering integral 


J2 j dpI{f,F} = 0 


( 21 ) 


0 ) 


we get for the drag current Jd (cf. with Ref. 

J, = -2eLV [ 

„ Jo 


dp 
2n h 




( 22 ) 


With the distribution function given by Eq. (TT2j) we have 


J d = -2eNVlL^[£- h j 


d 3 q 2n 


{I (P, n, Pe z \U\p + q z ,n, Pe z - q) | 2 (23) 


/ 0 2'rrh J (2 toK) 3 h 

X5(finp T Epe z £n,p+q z E Pe z — q) f npif fn,p+q z ) 

- |(p,n, Pe z + q\U\p + q z ,n,Pe z )\ 2 S(e np + E Pe z+q £n,p+q z Ep Sz )f n p+q z (1 fn,p)} 

where e z is the unit vector along the z axis. In the first term under the integral we change 
q —> — q and shift the integration variable p by q z , then 

d 3 q f qz dp 2n 


--^S 4 ks: 


27 rh h 


| (p,n,Pe z + q\U\p + q z1 n, Pe z )\‘ 


(24) 


so that the drag current is 
2MS R 


Jd — —Jo 


mn 2 K 2 


X5(e n p T E/p e ,_j_q E U) p—q z Ep Bz ) f np —q z (1 fn,p)i 


^ ./ d(lz io dp ^ np ~ q ^ l ~ ^ f ^^(*^> 19 * 1 ) 

x5 \_q\ — q 2 z (M/m — 1) + 2 q z (Mp/m + P)] . 


e(2eei) 2 LNmS-R 
(1 + k) 2 tt h 3 


Here we introduced 
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and a dimensionless quantity g(q_i_, \q z \) according to 


Sr9( q±> Wz\) = 


dRj_drj_e 


—ici±Hj_/h 


\^ n (r±)\ 2 K 0 (\qz\\^r ± \/h) 


where Sr is the cross section area of the ion beam. 
We get 

2MSr 


Jd ~ Jo 


[ dq z [ dpf np - qz (l - f nj) ) f dq ± g(q ± ,q z ) 
m7T n ~ n do do d 

x [?1 — q 2 z {M/m — 1) + 2 q z (Mp/m — P)] 

S [q]_ - q 2 z (M/m - 1) + 2 q~(Mp/m + P)] } . 


(25) 


(26) 


A. Linear response 

In the linear response regime 

V < T/p nF where p nF = \j2m\p - e r 


(27) 


the difference of (5-functions in Eq. (1261) can be expanded as (we again take into account that 
M/m 1) 


m 


2M\q 


<5 [q F — q 2 M/m + 2 q z [Mp/m + P)] — S [q F — q 2 z M/m + 2 q z (Mp/m — P)] 
{<5 [ mq]_/2Mq z — q z /2 + p + mV] — <5 \mq\/2Mq z — q z /2 + p — mb] } 

777 / d 

= H HI /2Mq z - q z /2 + p\ . 


M\q z \ dp 


Then integration by parts gives 




00 rln W* 


dq_ 

Qz Jo 


dp / dq ± g(q ± ,q z ) 


d 


Using 


we have 


d_ 

dp 


xS [mq\/2Mq z - q z /2+p\ —f np -q z { 1 - f n , P )• 


fnp—q z {1 fn,p ) (1 fn,p)d{q z p PnF) T fnp—q z d{j> PnF) 


(28) 


(29) 


Jd = Jo-^mV ^ f dq z J dq±g(q ± ,q z ) 

x [<& ~ ZPnFqz + rragl/M] (1 - f n , qz - PnF ) + 5 [ql - 2p nF q z - mq\/M\ fn, PnF -q z } ■ (30) 
Eliminating the (5-functions we get 


2$ _ r f 

Jd = Jo-^rnV E dq± < g{q±,PnF + pi) 


(i ~fn ,Pl) , „ \fn,p 2 


Pi 


+ d(q-L,PnF +P 2 ) J 


P 2 
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where p± = \J p 2 nF — mq\/M and P 2 = \J P 2 nF + mq\/M. The expression for J d can be 
simplified as follows 


Jd J 


45, 


R 


ir 2 h 2 


mV —— [ 

^ PnF J 


dq± 


g(q ± ,2p nF ) 

e q 2 j2MT + x ' 


If the ion beam cross section is of a circular form with radius a we have 

2 


g(q F , 2p nF ) = 


2hJi(aq±/h)\ r , 2 


aq _l 




K 0 (2p nF d/h), 


(31) 


(32) 


where J\(x) is the Bessel function of the first order and d is the distance between the ion 
flux and the wire. 

Below we will discuss in more detail a special case where g( qj_, 2 p nF ) does not depend on 
q^. For instance, this is the case provided 


V MT C h/a. 


(33) 


Then for J\f = 1 we get 

J d = Jo 4 ln ff— ( MT)—Kl(2p F d/K ), 
n z v F 


(34) 


where v F — p F /m is the Fermi velocity and 

e(2ee I ) 2 LNma 2 
0 = (1 + k) 2 H 3 


For the opposite case where 

VMT h/a 

the drag does not depend on temperature. For the values M = 10~ 22 g (Ga), T 
a = 10~ 5 cm this inequality can be easily satisfied. Then we have 


(35) 


(36) 

4K, 


Jd 


(¥) ' 


It is interesting to calculate the ratio Jd/ Ji for this case 

Jd e 32 (eej) 2 Lm 2 f 2p F d\ 
Ji ei (1 + hi) 2 7ih 3 vp 0 \ h ) 


(37) 


(38) 


Here one can use for Ko(s) equations (jHI) and (E]). 

For an estimate we assume the following values L = 10~ 4 cm, m = 7 • 10 _29 g, v F = 
2 • 10'cm/s, K = 10, p F d/h = 2, so that Kl(2p F d/h) = 1.3 • 10~ 4 . Then for J/ = 10~ 8 A one 
gets Jd = 2 • 10~ 9 A and the corresponding drag voltage Vd is about 


V d = 20 /jV. 

Naturally, if J/ goes up V d also goes up in proportion to J F 


(39) 














B. Nonlinear case 


We consider the simplest case of low temperatures assuming that 

V" > T/p F . 


(40) 


In our further calculation we will assume T — 0, then the integration due to the Fermi 
functions in Eq. (12611 is restricted and we get (the first or the second 5-function contributes 
for V > 0 and V < 0 respectively, so that the drag current changes its sign with V as it 
should) 

1 2 *—/ i i i,i i ~ i i i/I/,, i 

dp 




' PnF 


Qz 


%+r ^ 


' PnF 


* 2 PnF ( 1 Z d q z —p nF 


x j dqxg(qx,qz)6 [p - mV - q z (l - m/M)/2 + mq\/2q z M] 
The result valid for 


is 


Jd J 0 


ttH 2 


V < vp/2 

Y [ dq ± G[4Pp nF - q]_] f — g(q±,q z ) 

„ J Jv~ Qz 


(41) 

(42) 


where 0 is the step function and 

P± = PnF ± mV + \J ( PnF ± mV) 2 =F mq-jM 

(other cases are considered in Appendix A). 

For V <C v nF the integration variable q z is in the vicinity of 2p nF and we have 


a 2 f 

Jd = J °2^^ J d ^M^PPnF 




mV m q\ 


p nF M p 2 nF 


g(q±,2p nF ). 


(43) 


Eq. (l43l) substantially simplifies provided g(q_i_, 2 p nF ) does not depend on qj_; this is the case 
provided the ion flux cross section characteristic width a obeys the inequality 


Then 


where 


yjPp nF a/h -C 1. 

(44) 

Jd = J\ ^ 1 

n 

(45) 

(2 mV) 2 a 2 M 
' 1 “ 0 2 h 2 m ' 

(46) 


This expression is valid for V > 0, i.e. when the ion flux is directed ”to the right”. 
Then the momentum transferred to the electron system in the wire is also directed to the 
right and the current (since e < 0 ) flows in the opposite direction regardless of the sign of 
dragging ion charge. 
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Assuming that the distance d between the ion flux and the wire is much bigger than the 
characteristic cross section length of the wire and the flux we can write 


.g(q ± , 2p nF ) ~ Kl(2p nF d/h ), 


and 


Jd J1 E Kl{2p nF d/h). 

n 

Using for the function K 0 the approximate equation (J9]) we get 

J d = J iV] T—j e - 4knFd , 

knpd 

n 

where k nF = p nF /h. 

In the case ay/Pp nF /h 1 the drag current is linear in V 

J d = 4 J 0 MV Y, —K 2 0 (2p nF d/h). 

PnF 


(47) 

(48) 

(49) 


(50) 



FIG. 2: Drag current dependence on the velocity of the ion beam. We take Vj? = kv^ and 
k = 3/2. The first peak corresponds to V/v^ = 1/2 k (i.e. V/vp' 1 = 1/2) and the second peak 
corresponds to Vjv'p = 1/2. Here Jo = e(2ee/) 2 LNma 2 /(1 + K) 2 h 3 . 

III. CONCLUSIVE REMARKS 

We have developed a theory of Coulomb drag of electrons in ID ballistic nanostructure 
by an ion beam. This provides an example of drag of quasiparticles of the nanostructure 
by particles of the beam. It is worthwhile to mention that such a beam may consist not 
only of heavy ions but also of electrons. The free electron mass is usually bigger than the 
effective mass of conduction electrons, so that the adopted approximations of our calculation, 
M 77i, may remain valid in this case too. 
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The experimental setup should permit one to vary the velocity V within rather wide 
limits. We see however that to achieve a large drag effect one should choose the value of 
V near to v Fn /2 (see Fig. 2). This means in particular that the ion beam drag may be a 
useful instrument for nanostructure spectroscopy: it may make it possible to measure with 
appropriate accuracy the Fermi velocities v Fn in each channel n. 


Appendix A: Evaluation of the drag current for various ratios of a = V/v F 


We introduce dimensionless parameters 

a = mV/p nF = V/v nF 

and 

b = mqi/p 2 nF M 

and write q instead of g(l — m/M ) in the argument of the 5-function 

q poo pq 

dp + dq dp 

J 2 Jq-1 




X- / dqj_g{ja nF qx,p nF q)S [p - a - q/2 + b/2q] 

d J 


where 


Jq = e(2ee/) 2 LNma 2 / (1 + n) 2 h 3 . 

For a < 1/2 we get 

■h = Jo Y, (^1~) ^ J dq ± @[4a - b\ 


f A+(a,b) dq 9(PnFq±,PnFq ) 
A_|_(—a,— b) 0. 


where A±(a, b) = 1 + a ± y/(l + a) 2 — b. 
If 1/2 < a < 1 we get 


Jd — Jo 


„ V n J 77 J [ JA+(a-l,b 


pA + {a,b) 

+@[6 — 2a + l]0[4a — b] / dq 

J A+(—a,— b) 


dq 

- 1 , 6 ) 
g{PnFCj±,PnFq) 

q 


(Al) 


(A2) 


(A3) 


We will not give here explicit expressions for larger values of V/v F but rather present the 
simple expression for the drag current valid for h/a •C mVy/M/m 

Kl(2p F zd/h ) f 1 1 


Jd — 4Jq 


dz- 


J\z-a\ 2 -l)p 2 F /2mT _|_ ^ e ([z+a] 2 -l)p 2 F /2mT _)_ ]_ J 


(A4) 


Ja z \e K 

This expression is reduced to Eq. (1371) and Eq. (j5Uj) in the corresponding limiting cases. For 
mV T/v f the difference of the Fermi functions restricts the integration region so that 
for a < 1/2 we have 


Jd=iJo r d jM^iR 

Jl-a 


z 


(AS) 













1 + 0 ! 


Kl{2p F zd/K) 


(A 6 ) 


■Id — 4 Jo 


dz- 


for a > 1/2. The drag current calculated according to these simple formulas practically 
coincides with that calculated from the exact expressions and presented in Fig. 2 for the 
case J\f = 2 . 


Appendix B: Preferred velocity of the beam 

Let us differentiate Eq. (l26l) with respect to mV and determine the sign of derivative. We 
get under the sign of integral the following sum of 5-functions 


(Bl) 


(B2) 


— ^ {(5 [q 2 ± m/M - q 2 + 2 q z (p - mV)] + 5 [< q]_m/M - q 2 z + 2 q z (p + mV)] } . 

We integrate over p by parts and get (we denote the ratio J d 7r 2 h 2 / (2J 0 Snm) by j) 
dj f f qz 

^y = l dqz J rfq±£/( ' q± ’ Qz ^ J dp { 6 W± m / M ~ - ™V)] 

+ 5 [q\m/M -q 2 z + 2q z (p + mV)] } [f np -q z { 1 - fn, P )] 

/ OO p 

dq z J dq ± g(q ±1 q z ){ 1 - f q J {5 [q\m/M + q 2 - 2q z mV] 

+ S [q 2 ± m/M + q 2 z + 2 q z mV] } 

We take into account the strong Fermi degeneracy of the electron system so that 1 — f 0 = 
0, f 0 = 1. Using Eq. ([29j) we get 

dq±g(q±, q z ) {<5 [q 2 ± m/M + q 2 - 2q z (p nF + mV)] + 

5 [q 2 ± m/VI + q 2 z - 2q z (p nF - mV)] } 

poo p 

- dq z dq ± g( q ± , q z )5 [q 2 ± m/M + q 2 z - 2q z mV] 


dL = r 

w K„ r 


(B3) 


'PnF 


We again use fq z - PnF = 1 for q z > p nF and take into account that V > 0 so that the last 
5-function in Eq. f1B2li does not contribute. The second 5 -function in the first integral in the 
previous expression does not contribute as well and we arrive at 

dj 


dV 


= I+- 


where we introduce notations for the positive and negative integrals 
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If V/v n p < 1/2 we get 


U = 




g(q±,Qi) 


ql<AM P nFv 2 y/( p nF + mV) 1 2 3 - mq\/M 


(B6) 


where q\ = p n p + mV + a/( p n F + mV) 2 — niqjjM and /_ = 0 and the drag current is an 
increasing function of the beam velocity V. 

The integral /_ has nonzero values only if Vjv n p >1/2. If Vjv n p < 1 we have for this 
integral 


I 


dq± 


g{q±,Q2) 


q 2 ± <Mp nF (2mV-p nF )/m 2^(mV) 2 - mq 2 JM 


(B7) 


where = rnV + (■ mV ) 2 — mq\/M. /_ in this region becomes larger than I + (the latter 
being practically zero due to exponential dependence on q\) and the drag current turns into 
decreasing function of the velocity V. 

If V/v nF > 1 


/_ = 


q\<Mp nF v 


' q\>Mp nF (mV-p nF )/m 


dq ± 

dqp 


g(q±v-2) 

2y /, (mV) 2 — mq]_/M 

^(q± ; g3) 

2^(mV) 2 - mq 2 jM 


(B8) 

(B9) 


where q^ = mV — \J(mV) 2 — mq\/M. 

Therefore we see that the drag current has a maximum as a function of the beam velocity 
in the vicinity of V = v n p/2. 
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